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Abstract 



An effective black- hole-like horizon occurs, for electromagnetic waves in 
matter, at a surface of singular electric and magnetic permeabilities. In a 
physical dispersive medium this horizon disappears for wave numbers with 
k > k c . Nevertheless, it is shown that Hawking radiation is still emitted if 
free field modes with k > k c are in their ground state. 
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Hawking discovered that black-holes emit thermal radiation with temperature T# = 
Mp lanck /8iiM BH . For large black-holes, with mass M BH 3> M Planck , spacetime curvature 
is extremely small near the horizon, and the Hawking temperature is much smaller than 
the Planck temperature. It can be argued therefore, that as long as Mbh ^> Mpi anc k, 
the Hawking effect is essentially a low energy process, which is decoupled from quantum- 
gravity effects. Nevertheless, the standard derivation fl|] relys on free field theory on curved 
space-time, and describes the Hawking radiation as originating from vacuum fluctuations 
on scales exponentially smaller than 1/Mpi anc k On such scales, free field theory may 

not be a valid description, and quantum-gravity may become important. A question arises, 
therefore, to what extent is the Hawking radiation dependent on the details of the short 
distance structure of quantum gravity? 

In the absence of a theory of quantum-gravity, guiding toy models can be helpful. Unruh 
has suggested using a sonic analog of a black-hole 0. In his model, an effective event 
horizon is formed for sound waved on a flowing fluid background, on the surface where the 
flow becomes supersonic. Unruh has shown, that when a natural cutoff is implemented, sonic 
black- holes still gives rise to Hawking radiation In his approach the cutoff is implemented 
by modifying the dispersion relation for sound waves at high frequencies. This, in turn, alters 
the motion of modes with frequency close to the cutoff scale and gives rise to a new type of 
trajectories, which approach the horizon but eventually "reflect" back to infinity. It is not 
clear that a similar process is indeed realized for real black-holes. Further investigations, 
based on Unruh's model have been carried out in Refs. J^H . 

In order to gain further insight into the problem, in this letter, we provide and examine, 
a different toy model, which is based on phenomenological electrodynamics in matter. It is 
known that Maxwell's equations in a medium are formally analogous to electromagnetism 
in curved space-time |TIJ. In this case, an effective event horizon for electromagnetic waves 
in matter occurs on a surface of singular electric and magnetic permeabilities. As we shall 
show, formally, thermal radiation can still be emitted, if the initial state of the field is chosen 
appropriately. Nevertheless, this requires, non-physical, singular electric polarization on 
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arbitrary short scales, far beyond the scale where the electromagnetic field and the medium 
still couple. The short distance modifications that must be introduced in this toy model 
are clear: at high momentum or frequency the medium becomes dispersive and gradually 
the electromagnetic field and the medium decouple. This implies that for sufficiently short 
wave-lengths, the effective event horizon disappears. We shall show however that inclusion of 
such dispersive behavior does not eliminate Hawking radiation. The process which generates 
Hawking's radiation in the present toy model, differs however from that found for the sonic 
model. In the following we adopt units such that % = Kb — c = 1. 

To begin with, let us consider the macroscopic Maxwell's equations in a medium: 

V-B = VxE+f = (1) 
V-D = VxH-f = (2) 

where D and H are determined by the local electric and magnetic polarizations: D = E + P 
and H = B — M. We shall consider the special case of linear constitution relation of the 
particular form: 

H(r, t) = azB{r, t), D(r, t) = (az)^E(r, t), (3) 

i.e., the electric and magnetic permeabilities are e(z) = fi(z) = l/(az). 

With the gauge choice <p = and V ■ eA = 0, we obtain for the vector potential, the 
equation of motion 

.A. 

a 2 zV xzVxA + — - = 0. (4) 
ot 2 

Thus, although space-time is flat, with respect to the rest frame of the medium, waves 
propagate in an effective Rindler geometry: 

^effective = o? z 2 dt 1 — dz 2 — dx 2 — dy 2 . (5) 

The latter describes the space-time geometry near a black-hole horizon, of mass Mbh = 
(4a) _1 , for r /2Mbh *C 1 [11|. The event horizon corresponds here to the singular surface, 



z = 0, where e and /x — > oo. Although this model clearly breaks-down near the horizon, 
let us first examine its effect on electromagnetic waves and defer the discussion of a more 
physical realization to the sequel. 

The Hamiltonian of the system can be written as 

H = H I -H II = J™e~ 1 {z)(n 2 + {VxA) 2 yz-J \e(z) T 1 (d 2 + (V x A) 2 ) dz, (6) 

where A and D are conjugate canonical coordinates. The energy of a wave is thus negative 
in the region z < 0. This indeed suggests, that like a black-hole, the present model can emit 
thermal radiation. To see this let us quantize the field. 

For simplicity, we shall consider in the following only waves in the 5-direction, which 
correspond to radial transverse waves in the black-hole case. (Since we consider only the 
case t/2Mbh « 1, the vector nature of electromagnetic waves is unimportant.) The 
solutions of the wave equation (^) can be written as: 

9I (K)= ^ (7) 
[ 0, for z < J 

f 0, for z>0] 

= { } (8) 

I 1 giw^flnkl) f or z< oj 

where cj k = |ac|. The modes gi(n) correspond to an outgoing (right moving) wave for k > 
and ingoing wave for k < 0. Similarly, gni 1 ^) corresponds to an outgoing wave for positive 
k and an ingoing wave for negative k. The modes (0,|J) are orthogonal under the scalar 
product: 

/+oo 
t(z)g*i d t gidz, (9) 

and can be used to express the vector potential as 

/+oo 
(gi(K)ai(K) + gn(n)aii(n) + h.c.)d«. (10) 
-oo 

Here, ai and an are the annihilation operators in the domains z > and z < 0, respectively, 
and [a/-, an] = 0. 



Next, in order to determine what will be seen by a stationary observer at z > 0, we 
must specify an initial condition for the electromagnetic field at the effective horizon. Let us 
define the analogue Kruskal coordinates by U — —ze~ at and V = ze at . In terms of the U, V 
coordinates the effective metric is conformally flat and the solutions to the wave equation are 
outgoing right moving waves, exp(— iull), and ingoing left moving waves, exp(— iuV). The 



Unruh initial condition |12| amounts to the requirement that at the past horizon, V = 0, 
the outgoing modes, are in their ground state. Thus, time is defined with respect to light 
ray trajectories, z = —Ue at , with constant U. 

The analyticity of the outgoing modes, in the lower half of the U complex plane, implies 
that the combinations 

/!(«) = e™t** gi (K) + e-™l 2a g* u {-K), (11) 

and 

/ 2 («) = e™/**g n ( K ) + e^^gii-K), (12) 
have only positive frequencies with respect to [/-modes [O. It then readily follows that the 



outgoing state is thermal with (global) temperature 

r = s- (13) 

As in the standard derivation of Hawking radiation, this model involves ultra-high mo- 
menta. To see this, consider an outgoing photon described by the wave packet 

As we evolve the wave backward in time, the wave packet approaches the horizon like 
z ~ exp at, and gives rise to ultra-high momentum k w u exp(— at). Clearly this phenomena 
is also related to the occurrence of singular polarization on the horizon. 

Nevertheless, since the present model involves known physics, we do have some handle 
on these problems. When the wave length becomes comparable to the "molecular" length 



scale, we will have to account for the dispersive properties of the medium, which so far 
have been ignored. In particular, for sufficiently large wave numbers, the electromagnetic 
waves decouples from the medium degrees of freedom, and the effective geometry for the 
electromagnetic field becomes in this limit flat. Since in our model, the momentum k 
increases near the horizon, while the frequency, u, remains constant, we will introduce 
spatial- dispersion into the constitution equations @. As we shall see, this indeed eliminates 
the above mentioned difficulties. 

In order to describe both the inhomogeneity (in the z direction) and spatial dispersion, 
we shall now assume the following constitution relation for the electric field: 

D(z, t)= J e{z, z - z')E(z', t)dz . (15) 

and a similar relation, with [i(z,z — z') = e(z,z — z'), for the magnetic field. The first 
argument in e(z, z — z') allows for spatial non-homogeneity of the electric polarization, 
which in turn induces a low-momentum effective Rindler geometry. The second argument 
gives rise to a spatial non-locality around the point z, which is taken to be homogeneous, 
i.e. to depends only on z — z'. This allows us to simplify eq. (|T3p by going to the momentum 
representation. Assuming further that e is analytic in its first argument, we obtain: 

D(z,t)= J e ikz e(id k ,k)E(k,t)dk, (16) 

where 

e(id k ,k) = J e- ikz e(id k ,z)dz, (17) 

and E(k,t) is the Fourier transform of E(z,t). 

Using a similar relation between the magnetic field and the magnetic induction, we are 
led to the Hamiltonian 

H = J (^D(k,t)e-\id k ,k)D(k,t) + kA(k,t)e-\id k ,k)kA(k,t)^dk. (18) 

Here, D(k,t) and A(k,t) are the canonical coordinates. The equation of motion for the 
vector potential becomes 
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e-\id k ,k)k£~ 1 (id k ,k)kA(k)-u 2 A(k) = 0, (19) 

where A(k,t) = exp(±iut)A(k). Eq. (|I9|) can be further simplified, by noting that the 
solutions satisfy 

e- 1 (id k ,k)kA = ±uA. (20) 

The latter has a natural interpretation of an "eigenvalue-like" equation for the differential 
dielectric operator e~ l . 

We shall now make some assumptions on the behavior of e^ 1 . For k << k c , where 
l/k c is the "molecular" length scale, we must recover the non- homogeneous dependence 
e = 1/az, thus: e~ 1 (idk, k « k c ) —> iadk- On the other hand, for k >> k c , we expect the 
electromagnetic field to decouple from the medium, thus: e^iidk, k » k c ) — > (e c ) _1 , where 
e c is a constant. Assuming that e -1 interpolates smoothly between these two limiting cases, 
we shall represent the dielectric operator by 

e-^idk, k) = 6(k c - k)iad k + 9{k - k c )e~\ (21) 

where 9{k) = ~[1 + tanh(/c/A)] is the Heaviside step-function, and A determines the width 
of the transition. 

With this choice the solutions of eq. ( |T9"| ) are given by: 

G v (k,t) = ^exp{-^ [ k \ e ^ k '" k ^ A - ^(1 + jW-W^dkf)^ 



lexpj — [e 2 (*-M/A _ e 2( fc0 -M/Al _ il^k/kA _ t l e -^/*\Ei( \ k) - Ei(|fc 
k 1 2ae c 1 J a a 1 A A 
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(22) 

where Ei(x) = f^ OQ (e t /t)dt is the integral-exponential function, and u) n = \r}\. 

For k « k c (|22| ) reduces to r exp(±iu; r; t) exp(— In Substituting this into the 
Hamiltonian, we find that in the limit of k « k c , rjk > and rjk < correspond to positive 
and negative energy modes, respectively. In the other limiting case, k >> k c , we have free 
waves with u = \k\/e c , and the energy is positive for uj > 0. 



In order to relate these solutions with the modes (|7||f), and in particular with the thermal 

modes (|11],[1^), it will be useful to define the following low momentum modes: 

f -J= jfe-ie-**/*-*?^* for k > 1 
Mrj) = , (23) 

[ 0, for k < J 

f 0, for fc > ) 

I _^ = ^- 1 e to "*-^s ln l fc l), for fc < I 

for positive 77. For negative rj we replace above u — ► — cu. These modes form in the limit 
fc c — > oo, an orthogonal set under the scalar product: 

/+oo <-> 
dkg{e{id k , k) d t g 2 . (25) 
-oo 

For finite k c , we need to replace gj jj with the exact solutions G v (k,t). In the domain 
\k\ < k c , however, we can now expand the field as: 

/+oo 
(gi(v)Mv) + 9n(rj)an{rj) + h.c.)efy, (26) 
-oo 

where aj and ajj are the annihilation operators for k > and k < 0, states respectively. 

Next let us notice that although the spectrum of gi or gu in eqs. (|7|||) contains both 
positive and negative momentum k, the Fourier transform of the positive U frequency com- 
bination, fi in eq. flTTD, is: 



f°° he' ikz dz = 2*sinh(^)r(^ + l)^^), (27) 
J-oo v a ' v a J 

where k = r] > 0, and T is the Gamma function. The Fourier transform of ji yields gu. 
Consequently, gi and gu correspond to positive frequency [/-modes. In the limit of k c — > oo, 
we conclude that outgoing radiation modes, fi and / 2 , correspond in momentum space, to 
the modes gi and gu. 

When dispersion in high momentum is present this correspondence is no longer exact. 
gi and gu will be modified near the horizon in the regime \z\ < l/k c . Nevertheless, let us 
consider an outgoing wave packet, and a time to such that the wave packet is localized in 
a domain sufficiently far away from the horizon. At t — to we can describe the outgoing 
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Hawking photon and its negative energy pair by the mode /i, up to a small correction. The 
latter can be estimated by noting that in the momentum representation, large momentum 
amplitudes are suppressed like 1/k. Thus, if we instead start at t — to with the mode 
G v (k) gi, with rj > 0, it will correspond to fi up to corrections of at most 0(X c /M B h) 
where A c = 2n/k c . Notice that unlike f\, which is composed of both positive and negative 
parts, the corresponding mode G n {k) has a positive energy, i.e., an outgoing pair of thermal 
radiation are described by G v (k) modes in their ground state. 

To proceed with, let us now represent an outgoing Hawking photon and its pair at t — to, 
by a wave packet of the exact modes G v (k) with 77, k > 0: 

A g (k,t) = -^= r e^-^ )2/52 e- iut G v {k)duj. (28) 
V27r Jo 

t is chosen such that A g (k,t ) is mostly peaked at k < k c , and u c is taken to be much 
larger then u . We can then follow A g (k, t) as it evolves backwards in time, and observe the 
effects of dispersion. As is shown in the following figures, the momentum gradually increases 
and the wave packet, from which the radiation must have originated, is sharply peaked with 
k ~ k c . 
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Figure 1. The imaginary part of kA g (k, t), at t = 0, with k c = 20, uo = a = S = 1, A = 
0.5, e c = k c /uJo 
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As can be seen from Figure 1, at this stage the dispersion causes only a minor modification 
at k k c = 20, where the amplitude drops down to zero. 
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Figure 2. The wave packet at t = —5. 
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Figure 3. The wave packet at t = —10. 

In the last figure, the wave packet is completely localized within ~ 3 standard deviations, 
away from k c . It corresponds therefore to a wave packet of free waves which interact very 

10 



weakly with the medium. As we evolve A g even further back in time, the wave continues very 
slowly to shift towards larger k, and its width, Ak, narrows. For example at t — —40, the 
wave is localized at k ~ 22 with Ak « 0.2. The corresponding spread in space is Az ~ 1/ Ak. 
Therefore, the wave packet which gives rise to Hawking radiation, is at sufficiently early 
time much larger than the size of the black hole. 

The value of e c in the simulations above has been chosen as e c = k c /ujQ. Since the 
transition occurs at k « k c , and since for k > k c the frequency, ujq = k/e c , is fixed by the 
Hawking temperature, this choice corresponds to a minimal mismatch between the maximal 
dielectric constant, 

^max ^ 1/ ^^min ^ k^J ol, for low momentum waves, and 6^ for k ^> k c . 
This choice implies that the velocity of light (or g 00 ) drops down to a minimal value, l/e c , 
and remains constant for higher momenta. For other values, say e c = 1, free waves at k ~ k c 
have frequency uo ~ k c which does not correspond to the given Hawking frequency uoq. in 
this case we find that the Hawking radiation originates from a wave packet in the same 
range of k and the same amplitude as above, but which is strongly oscillating as a function 
of k. Although only weakly interacting free modes with uo ~ k c contribute, in the domain 
k > k c , this still effectively gives rise to the exact mode G v with low frequency uo = uo . (A 
related phenomena which acts in the reverse direction, yielding an effective high frequency 
mode from low frequency modes, was described in l i~3|| . We shall further describe this case 
elsewhere). 

Before concluding let us comment on some open issues. First we note that although the 
inclusion of spatial dispersion has eliminated the singular electric polarization at z = 0, we 
still can have e ~ k c /a ~ X c /Mbh >> 1- In other words, the use of linear constitution 
relations is questionable. Another question concerns the role of temporal dispersion, which 
has been so far ignored. In the present model, the latter does not seem to be essential 
in order to generate the Hawking radiation, since uo « k c at any time. Nevertheless, if 
we do assume a cutoff scale lu c , then the dielectric constant will typically be e(t — t' , z) ~ 
(za)~ l sin(c<j c (t — t'))/uo c . Therefore, temporal dispersion induces a non-locality, t — t' ~ 
l/azu c , which greatly increases near the horizon. This phenomena may be relevant for 
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understanding the final stages of black hole evaporation. 

In conclusion, we have established that outgoing radiation modes, originate in this model 
from nearly free modes, G v (k), with k m k c , which are weakly coupled to the medium. The 
latter are analogous to the Kruskal modes, exp(— iuU), in the standard picture, but do not 
contain wave numbers with k ^> k c . As a consequence, the initial state which gives rise to 
Hawking radiation, most naturally turns out to be the ground state of nearly free modes 
with k ~ k c , in an effectively flat background geometry. 
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